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ABSTRACT 

Empirical  kinetic  parameters  for  phenol-formaldehyde  polycondensate  materials  were 
surveyed  for  general  application  to  Air  Force  interests.  The  ranges  of  n.  A,  E  for  the  empirical 
kinetic  \nodel 


-(dw/dt)  =  w11  A  exp  (-E/RT) 

—1  16 

where  n=  0-5,  A  =  2.0x10  to  1.7x10  l/min,  E  =  3.7-72kcal/mole,  where  w  =  normalized 
weight,  t  =  time,  R  =  gas  constant,  T  =  temperature.  The  mass  spectrometry  and  thermo¬ 
gravimetry  thermal  analysis  methods  used  both  constant  heating  rate  and  isothermal 
approaches. 

Reexamination  of  data  with  anew  computer  program  gave  reduced  variations: 
n  =  1.56±0.07,  A  =  60.S±41.4,  E  =  10.38±1.20.  The  scopes  of  the  comprehensive  survey  and 
descriptive  thermogravimetry  computer  analysis  were:  (22  constantrate  +  3  isothermal)  vs  6 
constant  rate  experiments  by  8  vs  3  investigators  using  (5  cured  +  7  vitreous  fiber  reinforced) 
vs  4  cured  polycondensates. 

Parameter  inconsistencies  for  the  survey  were  partially  associated  with  an  atypical 
polycondensate,  laminate  and  polycondensate  particle  or  laminate  powder  samples,  or  a  reduced 
pressure  atmosphere.  Calculated  parameters  apparently  dependent  upon  the  two  latter  aspects 
were,  hot  averaged  with  other  values.  Small  curve-fitting  and  other  errors  potentially  enhance 
further  disagreement  as  they  often  gave  large  parameter  variations. 
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SYMBOL 
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c 

C 

e 

E 


exp 

f 

AH 

k 

K 

L 

1R(  ) 
logiQ(  ) 
m 

n 

P 


R 


TS(45°) 


w 


w„ 


X 


DEFINITION 

empirical  kinetic  parameter* 

system  concentration 

solid  specific  heat  at  constant  pressure 

base  of  natural  logarithms 

empirical  kinetic  parameter 

exponential  function  to  the  base  e 

function 

normalized  heat  of  reaction 
empirical  kinetic  parameter* 
integral  parameter 
thickness 

logarithm  to  the  base  e 
logarithm  to  the  base  10 
series  term 

empirical  kinetic  parameter 
pressure 
gas  constant 
surface  area 

error  variance  of  estimate 
population  error  variance  estimate 
temperature,  degrees  absolute* 
tensile  strength,  45°  to  fabric  warp 
reactant  (residue-free)  weight  fraction 
aggregate  sample  weight  fraction 
reactant  (residue-free)  weight 


SYMBOL 

UNITS 


l/min 


kcal/gm/°K 

dimensionless 

kcal/gm-mole 

dimensionless 

kcal/gm 

l/min 

dimensionless 


dimensionless 

dimensionless 

atm 

kcal/ gm-mole/°K 

2 

cm 

°K 

percent 

dimensionless 

dimensionless 

gm 
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a 

S,£ 

e 

SUBSCRIPTS 

SYMBOL 

c 

i 
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I 

m 

p.q 

r,s 

r 

1,2 

SUPERSCRIPTS 

SYMBOL 

I 

u 

nT 


SYMBOL 
UNITS 

dimensionless 
2  . 

cm  /min 
dimensionless 
°C 

DEFINITION 

calculated 

initial 

iteration  value 
isothermal 

maximum  rate,  maximum  value 
limiting  value  (q  span) 
limiting  value  (s  span) 
residue 
singular  value 


SYMBOLS  (CONT’D) 

DEFINITION 
oxygen  mole  fraction 
thermal  diffusivity 
increment 
temperature 


DEFINITION 
minimum  n  error 

mean 
d(  )/dt 
d(  )/dT 
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SYMBOL 

SYMBOL 

DEFINING  EQUATION 

UNITS 

a 

logi()r(AR/ET) 
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a' 

-w_  w~nt  exp  (X  ) 

Tm  ^  '  nr 
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al 

-w  =  kwn  exp  -  a2/L)] 

dimensionless 

a2 

cm 

a3 

-w  =  k  exp  [-a3(l  -  w)j 
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b 
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e  =  a  n  +bn+c 

*»  O  0  0 
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F(  ) 

2  ,  2 

V  )/se(  > 

G 

L(aCE/  AHRL)n+1 

n+2 .  .  n+1 
cm  /min 

g 

£  <woc  -  wo)2 
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J 

6,-9 

p-t-i  p 

°C 

k 

6,-9 

q+l  q 

°C 

P(X) 

X"1e“X  -  f°°  X'VXdX 

J 

dimensionless 

X 

q 

(0  -  0)/ j  +  i 

°C 

r 

X2eX  p(X) 

dimensionless 

s 

(0r  -  6a) /k  +  1 

°C 

se(a> 

se(wo)/q 
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se(b> 

%(%)/!  (1/T-17T) 

°K 

2.  x 

Se(wo> 

g'/(q  -  2) 

dimensionless 

2,  x 

s  (wj 
p'  o' 

[l(wo-^?2]/(q“1> 

dimensionless 

x 


AFML-TR-65-246 
Part  II 

FUNCTIONAL  SYMBOLS  (CONT’D) 


SYMBOL 

DEFINING  EQUATION 

SYMBOL 

UNITS 

w 

x/x.  =  (W0  -  wr)/(wo.  -  wp 

dimensionless 

X 

E/RT 

dimensionless 

y 

log10(K/T2) 

z 
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SECTION  I 


INTRODUCTION 

Part  I  of  this  report  considered  procedural  variables  in  the  constant  heating  rate  thermo¬ 
gravimetry  of  phenol-formaldehyde  polycondensates  {Reference  17).  The  best  sample  was 
an  adequately  cured  and  dried  molding  powder  free  of  any  postcure  surface  product.  Powder 
sieve  distribution,  pressure  or  weight  effects  were  small  to  500°or  600° C  for  runs  in  helium 
purging  flow.  Procedural  interference  increased  thereafter  to  about  4  percent  residual  weight 
and  250°  near  950°.  Fortunately,  a  change  in  a  procedural  variable  magnitude  gave  nearly 
parallel  thermograms.  Therefore,  the  residue-free  thermograms  were  relatively  free  of 
procedural  biasing.  With  this  result  a  kinetic  model  based  on  a  residue-free  normalized 
weight  was  adopted.  Empirical  kinetic  parameters  were  calculated  for  several  materials 
using  two  new  and  complimentary  curve-fitting  methods. 

The  residue-free  thermograms  and  elemental  compositions  of  fivepos toured  polycondensate 
powders  were  nearly  identical  from  300  to  810°  (Reference  37).  The  empirical  kinetic  model 
proved  representative  for  these  experiments.  Therefore  it  seemed  reasonable  that  similar 
empirical  kinetic  parameters  would  be  found  for  other  polycondensates  and  fiber  reinforced 
polycondensates.  On  this  basis,  a  survey  was  made  of  previous  constant  heating  rate  and 
isothermal  mass  spectrometric  and  thermogravimetric  thermal  analysis  studies. 

The  phenol-formaldehyde  materials  survey  was  disappointing  in  that  the  empirical  kinetic 
parameters  were  inconsistent  (Table  I).  A  systematic  attempt  was  made  to  explain  the  large 
differences.  First  a  survey  was  made  of  available  curve-fitting  methods  for  possible  computer 
program  use.  Although  no  entirely  satisfactory  technique  was  found,  a  program  was  evolved 
for  easy  and  precise  curve-fitting  using  a  dual  method  approach  developed  in  Part  I  of  this 
report. 

General  inspection  of  computer  program  or  manual  calculations  for  selected  experiments 
revealed  significant  differences  exceeding  the  normalizing  ability  of  a  residue-free  basis. 
Troublesome  aspects  v/ere  reinforced  laminate  and  polycondensate  particle  or  laminate 
powder  samples,  an  oxidative  atmosphere,  reduced  pressure,  and  unusual  behavior  of  the 
sample.  The  data  for  some  isothermal  thermogravimetric  and  mass  spectrometric  studies 
were  unsatisfactory  for  manual  calculations  and  were  not  amendable  to  the  computer  program. 

Data  clearly  disagreeing  with  likely  empirical  kinetic  parameters  for  a  cured  polyconden¬ 
sate  powder  or  unmanageable  data  was  rejected  for  this  study.  A  computer  program  analysis 
was  then  conducted  for  the  remaining  cases.  This  analysis  involved  six  constant  rate  runs  by 
three  investigators  using  four  materials  (Table  IX).  A  markedly  reduced  spread  in  the 
empirical  kinetic  parameters  was  encouraging  evidence  that  the  simple  kinetic  model  was 
probably  valid  for  most  cured  polycondensate  powders  and  reasonable  experimental  proce¬ 
dures.  Several  cases  of  model  inconsistency  or  significant  procedural  aberrations  were  also 
delineated  as  to  probable  causes. 

Computer  parameters  were  also  found  representative  for  seven  other  materials  either  by 
further  calculations  or  residue-free  thermogram  analogy  (References  4  and  37),  Further 
accepting  these  experiments,  the  survey  and  computer  study  scopes  were  then  respectively 
(22  constant  rate  +  3  isothermal)  vs  13  constant  rate  experiments  by  8  vs  4  investigators  on 
(5  cured  +  7  reinforced)  vs  11  cured  polycondensate  powders. 


*A11  temperatures  in°C  unless  noted  as  being  °F,  °K,  or  °R, 


1 


AFML-TR-65-246 
Part  II 


SECTION  n 

EMPIRICAL  KINETIC  PARAMETERS 

A.  EMPIRICAL  MODEL. 

Isothermal  kinetic  relations  which  empirically  represent  many  physicochemical  processes 
were 


-w  =  kwn  (l) 

k  =  A  exp  {- X )  (2) 

For  a  constant  rate  of  heating  T 

-w  =  ( Awn/T)  exp(- X  )  (3) 

B.  EVALUATION  METHODS  SURVEY. 

Numerous  methods  were  available  for  extracting  empirical  kinetic  parameters  from  experi¬ 
mental  thermograms.  Table  III  summarized  a  few  examples  for  constant  heating  rate  thermo¬ 
gravimetry  with  emphasis  on  thermosetting  polymer  studies.  The  approaches  were  categorized 
into  first,  one  or  more  necessary  heating  rates  and  then  residual  weight,  weight  loss  rate,  or  a 
combination  of  rate  and  weight  as  the  major  temperature  (or  time)  dependent  variable.  For 
example,  a  rearrangement  of  Equations  1  and  2  gave  a  method  using  rate,  weight,  and  an 
implicit  value  of  n 

lo,,0t-«w'n)  =  -r-~-E|0-  [y)  +  lo0|oA  (4) 

where  E,  A  were  respectively  evaluated  from  the  slope  and  intercept  of  a  parametric  plot  of 

log10(-ww  ”)  versus  1/T.  Other  popular  methods  used  a  difference  or  ratio  form  of  the  above 

equation,  an  integrated  version  of  Equation  3  for  a  1/T  parametric  plot,  a  maximum  rate 
of  weight  loss  theory.  Frequently  n  could  not  be  evaluated  directly  thereby  requiring  special 
consideration. 

The  precision  and  usage  of  the  different  methods  varied  considerably.  Critical  aspects  have 
been  described  previously  by  other  authors  but  will  be  briefly  reviewed  here.  Multiple  tem¬ 
perature  point  calculation  methods  had  the  advantage  of  being  self-checking  through  visual 
inspection  of  a  straight-line  parametric  plot.  Unfortunately,  such  plots  were  often  sensitive  to 
small  errors  and  thereby  biased  linearity  judgment.  This  difficulty  prevailed  near  the  beginning 
and  end  of  the  thermogram.  Measuring  weight  loss  rate  directly  from  the  thermogram  was 
imprecise  and  tedious  for  rapid  or  slow  rates.  Aids  as  computer  curve-fitting,  electronic 
differentiation  apparatus,  optical  devices,  large  and  accurate  records  served  to  reduce  error 
and  time.  Integral  weight  methods,  not  requiring  a  rate  input,  had  an  advantage  although  it 
was  usually  necessary  to  determine  n  by  separate  means. 

Multiple  heating  rate  methods  examined  kinetics  over  a  wide  temperature  and  time  spectrum. 
On  the  other  hand,  they  required  runs  at  different  rates  of  heating.  In  addition,  the  calculations 
proved  time  consuming. 

The  survey  of  analysis  methods  gave  nine  references  for  coherent  evaluation  of  the  three 
rate  parameters  n,  A,  E.  One  method  required  both  rate  and  weight  (References  22,  29  and  53). 
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In  addition,  the  parametric  plots  were  often  difficult  to  interpret.  Three  other  methods  used 
either  a  computer  program  or  complex  manual  calculations  not  otherwise  available  (References 
21,  36  and  44).  The  two  remaining  methods  specified  several  rates  of  heating  with  a  computer 
program  being  either  necessary  or  desirable  (References  2,  23  and  24). 

For  our  work,  no  single  method  was  entirely  satisfactory  for  easy  and  precise  curve-fitting. 
Therefore,  three  complimentary  computer  programs  based  on  the  Part  I  theory  were  developed 
for  this  purpose.  A  maximum  rate  experimental  method  (MAXRAX)  used  data  at  the  point  of 
maximum  rate  of  weight  change.  This  program  yielded  n  and  a  first  estimate  for  A,  E.  An 
integral  method  program  (TRIM)  used  corresponding  values  of  weight  and  temperature  from 
the  thermogram.  A  n  value  was  transferred  from  MAXRAX.  This  requirement  was  avoided 
for  the  third  program  (PRIM)  where  a  best  apparent  value  of  n  was  calculated  to  correspond 
to  the  minimum  of  an  overall  error  term  parabolic  fit  as  a  function  of  n  (Reference  36).  Tv;o 
additional  PRIM  features  were  a  comparative  pseudo-variance  analysis  and  two  temperature 
spans  for  A,  E  evaluation. 

C.  PHENOL-FORMALDEHYDE  MATERIALS  SURVEY. 

Approximate  ranges  of  experimental  conditions  and  empirical  kinetic  parameters  for  the 
survey  were  summarized  in  Table  I  with  detailed  results  in  Table  II.  In  most  constant  rate 
work,  5-400  mgm  of  a  fine  powder  sample  was  run  in  a  porcelain  crucible  to  about  950°  at 
3° -10° /min.  Bulk  material  was  frequently  pulverized  under  liquid  nitrogen  to  reduce  friction 
heating.  A  purging  flow  of  inert  gas  was  generally  use  with  a  few  moderate  and  high  vacuum 
experiments  being  reported.  Available  thermograms  and  commercial  processing  suggested 
a  desirable  minimum  cure  or  postcure  temperature  of  about  125°  with  temperature  being  more 
important  than  time  for  thermogravimetry.*  Therefore,  the  maximum  processing  temperature 
was  indicated  as  a  crude  thermal  history  index. 

The  empirical  kinetic  parameters  for  constant  heating  rate  thermogravimetry  of  cured  and 
vitreous  fiber  reinforced  materials  overlapped.  The  higher  n,  A,  E  parameters  were  for 
reinforced  samples.  The  largest  set  of  values  were  designated  temperature-dependent  and 
varied  with  each  of  three  curve-fitting  methods. 

The  parameter  E  for  the  single  isothermal  thermogravimetry  source  for  a  cured  poly¬ 
condensate  was  comparable  to  an  isothermal  mass  spectrometric  analysis  result,  with  n  =  0 
in  both  cases.  Constant  rate  mass  spectrometric  thermal  analysis  E  values  were  temperature- 
dependent,  ranging  from  16  up  to  48  kcal/mole,  the  highest  reported  for  a  cured  polycondensate. 
The  parameter  n  was  zero;  unfortunately,  A  results  were  not  given  for  this  study. 

The  empirical  kinetic  parameters  for  the  survey  were  somewhat  mutually  compensating  for 
different  experiments.  For  a  given  n,  an  increase  in  A  gave  an  increase  in  E  for  ea  1.  investi^ 
gator  with  only  four  close  exceptions.  For  a  given  n  and  experimental  set,  an  increase  in  T 
was  consistent  with  an  increase  in  both  A,  E  with  only  one  minor  exception. 

For  two  sets  of  A,  E  and  with  constant  n,  T  identical  values  of  w  would  correspond  to 

(  riA|E2)/(r2A2E|)  =  ( T, /t, )  exp  [  (E(  -  E2  )/Rt]  (5) 

For  Tj  =  Tg ,  there  was  a  single  intersection  of  th^  two  thermograms  and  large  differences 
were  suggested  elsewhere.  The  low  precision  at  hand,  however,  let  sets  of  A,  E  and 


♦Unpublished  work  by  the  author. 
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especially  n,  A,  E  approximate  the  same  thermogram.  Thus  various  parameter  sets  were 
dependent  upon  curve-fitting  sensitivity.  The  same  effect  clearly  existed  for  different  heating 
rates.  This  mutual  compensation  tendency  dictated  a  new  and  precise  curve-fitting  technique  for 
the  reanalysis  of  previous  data. 

A  factual  heating  rate  or  temperature  dependence  of  A,  E  was  reported  by  a  number  of 
investigators.  This  was  believed  to  be  due  in  part  to  mutual  compensation  from  inaccurate 
curve-fitting  as  further  illustrated  in  a  subsequent  section  of  this  report.  In  a  formal  sense 
Equation  3  was  invalid  for  temperature  dependent  empirical  kinetic  parameters.  The  use  of 
this  relation  and  associate  curve-fitting  techniques  gave  pseudo-values  although  the  errors 
involved  were  probably  small.  Although  at  least  one  partial  solution  for  temperature  dependent 
values  of  A,  E  has  been  obtained,  no  suitable  evaluation  method  could  be  derived  in  detail 
(Reference  20).  Further  factual  dependence  of  the  empirical  kinetic  parameters  of  the  invariant 
model  on  such  procedural  aspects  as  ambient  pressure,  sample  form  and  weight,  and  so  forth 
was  a  major  consideration  in  the  reexamination  of  thermograms  using  a  new  computer  program. 
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SECTION  m 

EXPERIMENTAL  MATERIALS  AND  METHODS 

Table  IV  summarized  the  major  materials  for  this  report.  Molding  and  other  samples  of 
Monsanto  RI4009  commercial  phenol-formaldehyde  polycondensate  were  previously  examined 
(Reference  17).  The  remaining  materials  were  identical  to  those  of  either  isothermal  thermo¬ 
gravimetry  or  thermal  expansion  study  (References  12  and  17). 

A  polycondensate  casting  and  molding  were  pulverized  under  liquid  nitrogen  with  a  jewelers’ 
cutter  (moly-steel  10,000  rpm)  and  sieved  to  a  standard  particle  distribution.  A  reinforced 
laminate  powder  was  prepared  in  a  similar  manner  but  using  a  hole  saw  (moly-steel,  3/4  inch 
diameter,  780  rpm).  Integral  glass  fabric  reinforced  samples  about  one  inch  square  by 
1/8  inch  thick  were  used  for  isothermal  air  aging. 

The  modified  Aminco  Thermo-Grav  recording  thermobalance  and  experimental  methods  for 
constant  heating  rate  thermogravimetry  were  described  in  Part  I.  Laminate  plate  aging 
experiments  were  conducted  in  a  forced-air  convection  oven  from  300° -800°  F  (about  150  - 
430°  C)  to  1000  hours.  The  procedures  and  mechanical  property  evaluations  have  been  further 
detailed  or  referenced  elsewhere  (Reference  19). 
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usually  agreed  with  the  “best”  n  value.  Small  analytical  and  perhaps  experimental  errors  at 
the  ends  of  the  thermogram  severely  affected  A,  E  results  as  the  w,  K  values  were  especially 
sensitive  here.  This  was  further  illustrated  by  Figures  8  and  9. 

The  parametric  study  further  suggested  that  visual  inspection  of  data  point  deviations  from 
a  mean  line  for  a  parametric  plot  was  not  reliable  for  the  precise  evaluation  of  empirical 
kinetic  parameters.  The  third  computer  program  PRIM  was  designed  to  provide  other  means 
for  proving  reliability. 

For  PRIM,  six  incremental  values  of  n  were  chosen  around  the  MAX  RAX  n  result.  Essential¬ 
ly,  TRIM  was  then  repeated  for  these  six  values.  An  overall  error  term  g  was  calculated  for 
each  case 

0  '  I{«oc_*o)2  02) 

A  best  value  n  ’  was  then  computed  as  a  minimum  for  an  approximation  to  g  dependence  upon 
n  using  a  parabolic  method  of  least  squares 

g  =  oQn2  +  bQn  +  cQ  (13) 

Then  the  result  n  ’  was  used  for  calculating  corresponding  values  of  A,  E  using  TRIM. 

PRIM  used  two  temperature  spans.  The  first  or  q  span  was  primarily  for  securing  the 
best  value  of  n.  The  second  or  s  span  was  used  for  calculation  of  A,  E  for  the  entire  thermo¬ 
gram.  The^  s  range  was  a  compromise  between  available  data  and  likely  pyrolysis,  140’  - 
1010°  in  30°  increments.  The  q  span,  which  also  used  3 (f  increments,  was  roughly  chosen  for 
each  case  using  results  from  an  initial  parametric  study 

(•op.  V  *  {0  98.°52)  04) 

(wpfwq)  *  (0.95,0.05)  (15) 

(0p ,  0q )  -  (320*880*)  (16) 

A  feature  of  PRIM  was  evaluation  of  four  indices  of  conformance  for  the  calculated  and 

experimental  thermograms.  One  was  g'  and  a  second  was  ±(w  -  w_)  ,  this  value  times 

oc  Or  m 

100  giving  the  maximum  sample  weight  percent  error  for  the  s  temperature  span.  Comparative 
2  2 

pseudo- vari ances  se( a),  se(b)  were  calculated  assuming  that  Equation  11  was  valid  for  a 

least  mean  squares  analysis  (Reference  42).  These  parameters  were  respectively  indicative 
of  the  deviations  of  A,  E  and  E  from  a  mean  (Appendix  II). 

Absolute  methods  were  not  available  for  either  the  linear  or  parabolic  approximating  least 
mean  squares  cases.  As  the  pseudo-methods  could  fail,  several  checks  were  made.  A 
coefficient  A’  was  calculated  by  an  independent  relation  to  confirm  the  value  of  A  from  Equation 
10.  Several  conformance  indices  were  also  calculated  for  the  parabolic  pseudo-method  as 
further  discussed  in  Appendix  II. 

The  programs  were  written  in  Fortran  IV  for  an  IBM  7044/7094  system.  Depending  upon 
the  program,  from  2-10  punched  input  cards  were  necessary.  Only  descriptive  outlines  were 
given  herein  as  several  sets  of  symbols  were  used  and  the  programs  were  not  optimized  for 
minimum  running  time. 
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SECTION  IV 

EMPIRICAL  KINETICS  COMPUTER  PROGRAMS 


A.  THEORY. 

The  general  theory  for  the  computer  programs,  developed  in  Part  I,  was  also  reviewed 
herein  as  Appendix  I  for  reference  convenience.  Sequential  steps  for  the  three  programs  and 
a  thermogram  construction  subroutine  were  summarized  using  flow  sheets  (Figures  1-4) 
and  outlines  (Appendix  II). 

The  maximum  rate  experimental  program  MAXRAX  essentially  solved  four  equations 
using  an  iterative  procedure  and  two  tables  * 


^ In  r (  n  *  —  1)  +  ij 

n  =  - - — ■  +  1 

ln  wm 

(6) 

A  =  w^“n{ix/nT)mexp(Xm) 

(7) 

E  8  -WTmR(nT2/")" 

(8) 

r  =  X*«X|”p<Xj 

(9) 

The  necessary  input  data  were  evaluated  at  the  thermogram  maximum  rate  of  weight  loss. 
This  point  and  the  associate  heating,  and  weight  loss  rates  were  visually  estimated  using  a 
straight  edge  and  an  optical  device  based  on  the  alignment  of  two  partially  reflected  images 
of  the  curve. 

The  second  computer  program  TRIM  solved 

=  Tririolr] +  **»'««'*"  (l0) 

which  was  derived  from  the  integrated  solution  of  the  constant  rate  kinetic  equation.  Other 
than  experimental  data,  the  inputs  n,  r  were  from  MAXRAX  and  the  approximation  T  =  Tm 
was  made. 

For  TRIM,  the  parameters  A,  E  were  calculated  for  an  ideal  equation 

y  -  bz  +  a  (II) 

by  the  standard  linear  method  of  least  squares.  This  was  an  approximating  pseudo-method 
as  all  term3  were  interdependent,  experimental  precision  was  variable  and  the  term  r  slowly 
decreased  with  increasing  temperature  (about  ±2  percent  maximum). 

Thermograms  for  -325  sieve  molding  powder  at  three  heating  rates  were  fairly  smooth  and 
consistent  excluding  the  final  three  percent  of  weight  loss  above  70 (f  (Figure  5).  A  slight 
discontinuity  near  400°  for  the  9?  /min  run  was  also  ignored.  Using  TRIM,  systematic  variations 
of  n,  A,  E  and  the  temperature  span  were  used  to  study  these  thermograms.  The  “best” 
values  of  A,  E  for  minimum  error  between  calculated  and  experimental  thermograms 
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separation  by  a  constant  sample  weight  difference  with  a  small  error  occurring  at  intersections. 
Thus,  inspection  of  Figure  7  revealed  a  residue-free  resemblance  to  about  500°  with  similarity 
for  powder  samples  to  900°.  A  particle  was  an  unsatisfactory  geometry  in  agreement  with 
other  results  (Reference  18).  Reference  10  further  notes  that  complete  drying  would  probably 
have  given  better  agreement  to  500°;  smaller  errors  to  900°  would  also  result.  The  low 
temperature  weight  losses  also  suggested  too  low  a  processing  temperature  with  some  solvent 
retention.  The  small  thermogram  fluctuations  were  abnormal  and  could  have  been  due  to  an 
instrument  artifact,  replotting,  or  other  unknown  factors. 

PRIM  empirical  kinetic  parameters  for  polycondensate  casting  and  glass  fabric  reinforced 
laminate  powders  were  not  the  same  (Table  VII).  The  laminate  powder  sample  was  more  stable 
for  a  residue-free  basis.  This  apparent  stability  was  believed  primarily  due  to  glass  particle 
sintering  at  the  higher  temperatures  with  encapsulation  of  otherwise  decomposable  material. 
For  example,  a  glassy  residue  with  brownish  inclusions  was  found  when  the  laminate  powder 
helium  residue  was  rerun  in  air.  Differences  in  resin  contents  calculated  for  the  casting  and 
laminate  powders  in  helium  (19.4%),  laminate  powder  in  air  (24.3%),  and  a  laminate  plate 
(27.0%)  were  not  unreasonable  in  that  any  encapsulated  material  would  be  ultimately  weighed 
as  “glass”  and  encapsulation  would  piobably  be  least  effective  in  air.  A  plate  resin  content 
of  27.0  percent  was  a  moisture-free  normalized  value  for  long  time  burnout  in  air  at  540° 
(Reference  6).  Previous  results  with  an  uncured  resin  and  heat  cleaned  18I-“E”  glass  powder 
(50/50  mixture,  -400  sieve)  also  suggested  encapsulation  at  temperature  above  about  500° 
(Reference  18).  The  helium  residues  near  950°  were  59.5  percent  for  the  mixture  versus  56 
percent  for  the  polymer  alone. 

It  proved  necessary  to  “dry”  the  casting  and  laminate  powders  prior  to  thermogravimetry 
(Table  VII).  The  conditions  were  120°  for  25  minutes,  which  gave  an  average  stabilized  weight 
loss  of  1.45  percent.  About  25  minutes  were  used  for  heating  to  120°. 

Empirical  kinetic  parameters  were  calculated  for  phenol-formaldehyde  polycondensate 
powder  and  sized  glass  fabric  reinforced  laminate  samples  (Table  VIII).  The  powder  sample 
was  more  stable  than  the  plate  for  both  helium  purging  flow  and  a  moderate  vacuum.  The  spread 
in  A,  E  was  fairly  small  for  all  four  experiments  with  n  values  being  similar  for  plate  samples 
but  depending  upon  environment  for  powders.  Enhanced  stability  was  found  under  vacuum  for 
both  materials.  This  has  been  attributed  to  an  increase  in  pyrolysis  product  fragmentation  with 
increasing  pressure  and  may  be  associated  with  inefficient  removal  of  pyrolysis  gases  by  a 
purging  flow  (Reference  18). 

PRIM  and  Reference  53  parameters  were  not  the  same  for  a  powder  and  plate  under  vacuum 
(Table  VIII).  The  thermogram  error  index  ±(wqc -wQ)m  suggested  that  PRIM  was  more  repre¬ 
sentative  for  the  powder  with  PRIM  superiority  being  conclusive  for  the  plate.  There  were 
consistent  parameter  variations  for  a  polycondensate  powder  versus  a  plate  geometry  as 
compared  to  a  powder  versus 'a  laminate  powder  (Tables  VII  and  VJII).  Plate  pyrolysis  has 
been  suggested  to  be  influenced  by  catalysis,  geometry,  glass  fabric  finish,  or  possibly  other 
effects  (Reference  53).  Polycondensate  and  laminate  powder  differences  were  believed  primai- 
ily  due  to  particulate  glass  sintering. 

Five  alkaline  catalyzed  polycondensate  powders  have  undergone  similar  gravimetric 
pyrolysis  once  curing  was  complete  (Reference  37).  PRIM  parameters  for  a  sample  coded 
B-C  were  comparable  with  other  results  (Table  IX).  The  necessary  thermogram  was  con¬ 
structed  from  tabulated  B  -  wq  values  and  extrapolated  to  950°  (Reference  37).  As  nearly 

identical  residue-free  thermograms  were  reported  for  these  materials,  the  sample  B-C 
parameters  were  also  representative  for  all  five.  They  included  three  catalyzed  by  Ca(OH)2 

(P/F  ratios  of  0.36,  0.42,  0.50)  and  the  commercial  products  CTL-91LD  and  SC1008, 
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PRIM  curve-fitting  indices  and  empirical  kinetic  p*  rameters  were  averaged  for  six 
experiments  using  an  inert  gas  purge  (Table  IX'-.  A  definitive  pseudo-variance  analysis  was 
intended  to  aid  the  selection  of  experiments  Lut  could  not  be  completed  for  this  report. 
Qualitative  results  follow: 

(a)  Molding  powder  (Tables  V  and  VI).  The  spread  in  MAX  RAX  and  TRIM  n.  A,  E  values 
for  three  heating  rates  exceeded  the  spread  of  the  reference  set.  For  PRIM,  the  3° /min 
parameters  were  adjudged  the  best  of  the  reference  set.  The  9° /min  results  gave  upper  A,  E 
limits  for  the  reference  group  although  also  giving  a  good  calculated-experimental  thermogram 
fit. 


(b)  Pure  polycondensate  and  glass  fabric  reinforced  powder  samples  (Table  VII).  The 
polycondensate  values  were  used  for  the  reference  group  subject  to  being  the  lowest  E  case. 
The  high  n.  A,  E  laminate  powder  results  were  questionable  even  in  view  of  an  excellent 
overall  fit  resulting  in  the  lowest  conformance  indices  of  any  run. 

(c)  Pure  polycondensate  powder  and  glass  fabric  reinforced  plate  samples  (Table  VIII). 
The  helium  purge  powder  parameters  were  adjudged  satisfactory  for  the  reference  set.  The 
vacuum  results  were  marginally  withheld  due  to  high  n,  low  A  values  as  well  as  a  poor  overall 
fit.  Although  similar,  the  plate  sample  values  for  both  helium  and  vacuum  were  not  acceptable 
because  of  a  low  n  and  poor  fit. 

(d)  Polycondensate  powder  (Table  IX) .  The  parameters  for  the  only  laboratory  synthesized 
material  were  used  in  the  reference  set  subject  to  being  the  maximum  overall  fit  error  case. 

(e)  Reference  set  materials  (Table  IX).  There  was  still  a  tendency  for  mutual  compen¬ 
sation  of  the  parameters,  especially  for  A,  E  between  experiments.  There  were  no  consistent 
trends  for  n.  A,  E  with  a  change  in  T.  There  was  a  tendency  toward  higher  values  of  g‘ , 

2  2 

se(a),  sg(b)(or  average  error)  but  not±(wQc  -  wo)m(ormaximum  error)  at  the  higher  values  of  T. 

Within  the  limited  spectrum  of  experimental  conditions,  no  consistent  dependence  of  the  con¬ 
formance  indices  or  parameters  upon  other  variables  as  grinding  method,  sample  weight,  and 
so  forth  was  evident.  Although  the  spread  in  the  parameters  was  low,  especially  when  compared 
to  the  initial  materials  survey,  fair  errors  resulted  if  one  thermogram  was  predicted  using 
either  the  reference  set  parameters  or  results  for  another  experiment.  Rough  estimates 
implied  ±(woc  -  wQ)m  of  about  0.03-0.06  over  a  narrow  temperature  range  with  up  to  0.03 

average  error.  Ihe  qualitative  selection  of  experiments  for  averaging  and  the  averaged  set 
was  biased  by  both  low  overall  accuracy  and  PRIM  limitations  to  the  degree  that  all  observa¬ 
tions  in  the  above  paragraphs  were  tentative  and  subject  to  future  revision. 

The  residue-free  weight  w  for  the  empirical  kinetic  model 

w  =  (wQ  -  wf  )/( wQj  -  wf)  (17) 

was  for  a  continuous  thermogram.  The  rate  -wqT  could  approximate  either  zero  at  two  points 

or  one  maximum,  but  could  not  be  a  minimum.  These  rules  permitted  quick  rejection  of 
noncomforming  thermograms.  For  multi-step  thermograms  (with  -w^  approaching  zero  at 

more  than  two  points)  the  initial  unity  value  of  w  .  would  have  become  w  for  each  successive 

oi  r 

step.  Equation  17  in  a  different  form  permitted  easily  constructed  conversion  diagrams  for 
graphical  calculations  or  rapid  estimates  (Figure  8).  General  effects  of  n,  w  upon  K  were 
often  pronounced  (Figure  9).  K  values  were  insensitive  to  n  but  strongly  dependent  upon  w  for 
w>0.8.  This  portion  of  the  thermogram  was  limited  for  n,  A,  E  evaluation  and  was  unfairly 
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biased  by  small  errors.  For  a  w  of  less  than  about  0.2,  K  was  strongly  dependent  upon  both 
n,  w.  Although  perhaps  attractive  for  n  evaluation,  small  errors  also  exerted  an  unfair  weight¬ 
ing  here. 

Accuracy  and  precision  levels  for  constant  heating  rate  thermogravimetry  were  reevaluated 
and  the  new  results  summarized  as  Table  X  (Reference  17).  In  addition,  slightly  revised  wf 
data  were  used  for  some  PRIM  calculations  (Table  V). 

MAXRAX  reflected  both  experimental  and  model  fitting  errors  as  primarily  compared  to 

an  overall  model  error  for  PRIM.  A  ±{w  -  w  )  difference  between  the  two  was  therefore 

'  oc  o'm 

roughly  the  experimental  weight  error  at  the  point  of  maximum  rate.  Calculated  values  were 
reasonable  and  varied  up  to  twice  the  estimated  experimental  precision  (Table  X).  An 
argument  similar  to  the  one  above  was  possible  to  show  that  the  difference  between  A,  A* 
approached  the  error  in  A  associated  with  the  experimental  error  at  the  point  of  maximum 
rate.  In  most  cases,  this  difference  was  encouragingly  small  (for  example,  see  Table  VI). 

The  program  and  past  studies  revealed  a  variety  of  troublesome  aspects  in  the  evaluation  of 
consistent  empirical  kinetic  parameters  for  phenol-formaldehyde  polycondensate  materials 
(References  17,  18  and  20).  Specific  examples  exceeding  the  normalizing  ability  of  residue-free 
empirical  kinetics  were: 

(a)  Critical  curve-fitting  aspects  for  a  single  thermogram  or  similar  set  including  ~ 

1.  Data  point  spacing,  reaction  rate,  or  temperature  range  incorporation  and  possible 
magnification  of  errors  (especially  near  thermogram  ends). 

2.  Improper  use  of  the  invariant  model  in  evaluating  temperature  dependent  parameters 
or  thermograms  with  more  than  one  maximum  rate  of  weight  loss,  and  so  forth. 

3.  Mutual  compensation  of  different  sets  of  parameters  permitting  approximate 
representation  of  a  single  thermogram. 

4.  The  very  flexible  nature  of  the  kinetic  model  for  approximate  representation  of 
different  thermograms. 

(b)  An  oxidative  or  reduced  pressure  atmosphere. 

(c)  Reinforced  laminate  or  polycondensate  particle,  laminate  powder,  or  polycondensate 
film  samples. 

(d)  Samples  contaminated  with  postcure  surface  product  or  solvent,  inadequately 
cured  or  dried,  or  too  large. 

(e)  Two,  three,  or  occasionally  many  stepped  thermograms.  This  atypical  behavior 
appeared  pronounced  for  large  sample  particles  and  perhaps  minute  thermobalance  oxygen 
contamination;  low  inert  gas  flow  rates  could  have  resulted  in  either  contamination  or  ineffi¬ 
cient  pyrolytic  gas  removal.  Stepping  was  perhaps  accentuated  to  varying  degrees  by  (b)-(d) 
items  as  well  as  the  lack  of  pressure  during  sample  polymerization. 
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SECTION  V 

TYPICAL  APPLICATIONS  OF  EMPIRICAL  KINETIC  PARAMETERS 

A.  EMPIRICISM  OF  n,  A,  E 

The  kinetic  parameters  were  strictly  empirical.  There  was  no  factual  significance  in 
terms  of  classical  physicochemical  mechanisms  (Table  XI).  Attempts  to  assign  primary 
governing  mechanisms  could  not  resolve  complex  intra-relations  between  the  sample, 
experimental  method,  and  possible  large  numbers  of  chemical  reactions  and  physical  processes 
during  pyrolysis.  The  kinetic  model  undoubtably  grossly  represented  a  variety  of  complex 
mechanisms. 

B.  RELATIVE  THERMAL  STABILITY 

The  empirical  kinetic  parameter  E  has  often  been  used  as  an  index  of  relative  thermal 
stability.  For  example,  this  parameter  was  good  for  many  thermoplastics  with  a  n  of  a  near 
integer  and  similar  large  A  values.  For  thermosets  with  small  values  of  A,  E  and  high  per¬ 
centage  of  final  residue,  the  three  parameters  n,  A,  E  proved  closely  interrelated.  More 
suitable  indices  were  the  maximum  rate  loss  rate  and  corresponding  temperature.  The 
normalized  weight  at  this  point  proved  largely  controlled  by  n,  r  at  about  0.5.  For  a  known 
n,  A,  E,  T,  the  maximum  rate  indices  were  only  available  by  interpolation  or  trial  and  error 
calculation 


II 

V6 

nw"  1  (AR/ET)  «xp(-X  ) 
m  m 

(18) 

Km  5 

<AE/RT)p(Xm) 

(19) 

~WTm“ 

w  (  X/nT) 
m  m 

(20) 

C.  HYPERTHERMAL  ABLATION 

Kinetic  models  similar  to  the  current  one  have  been  used  extensively  for  hypei'thermal 
ablation  calculations  and  computer  programs.  Three  regions  were  used  in  one  descriptive 
char-formation  ablation  analysis  (Reference  49) .  These  regions  were  a  carbonaceous  char 
layer,  a  reaction  zone,  and  the  virgin  material.  A  char-reaction  zone  interface  recession  rate 
for  the  interface  temperature  was* 


n+2  2n+l  -2(n  +  l)  -X 


(L)  =  AGn!  (  1  -  wr  )  X 

e 

(21) 

G  =  L(aCE/AHRL)n+l 

(22) 

For  a  given  G,  the  lowest  recession  rate  would  result  for  a  large  E,  wf  and  a  small  A,  T. 

Clearly,  small  uncertainties  in  the  parameters  could  give  a  large  L  error.  Some  typical  para¬ 
meter  values  and  results  from  another  computer  study  were  summarized  by  Table  XII  and 
Figure  11. 


♦This  derived  form  uses  the  reasonable  approximations:  1«AHRX;  T— 0  as 
L— oo  ;  (wall/pyrolysis)  mass  flux  « (1  -  wf).  Typically,  L  <<  1. 


12 


AFML-TR-65-246 
Part  II 


This  application  of  empirical  kinetic  parameters  to  a  complex  environment  appeared 
primarily  because  of  inadequate  experimental  methods.  The  simple  thermogravimetry  ex¬ 
periment  for  an  isolated,  nearly  isothermal  sample  was  extrapolated  to  an  internal  pyrolysis 
zone  subject  to  intense  and  rapidly  changing  catalytic,  compositional,  pressure,  stress, 
thermal,  volumetric  and  other  gradients.  For  example,  the  time  for  the  receding  steady  tem¬ 
perature  distribution  of  Figure  11  to  traverse  a  pyrolysis  zone  thickness  of  about  0.35  mm 
would  be  1.4  second,  corresponding  to  an“averageM  heating  rateofl8,500°/min.  In  other  cases, 
especially  at  low  heat  flux,  the  heating  and  recession  rate,  and  temperature  distribution  was 
less  abrupt  with  a  resultant  thick  pyrolysis  zone.  Excluding  limited  high  heating  rate  and  high 
temperature  work,  no  definitive  empirical  kinetic  model  or  parameters  could  be  found  for  this 
complex  environment  or  a  fair  experimental  simulation. 

D.  MECHANICAL  PROPERTIES. 

The  use  of  empirical  kinetic  models  and  parameters  for  the  analysis  and  correlation  of 
mechanical  properties  has  been  limited  (Reference  16).  In  one  study  with  this  intended  goal 
the  present  model  failed  to  represent  the  air  aging  weight  changes  of  CTL-91LD/Glass  Fabric 
(Figure  12,  Reference  19).  The  isothermal  thermograms  for  this  and  similar  materials 
suggested  two  gross  controlling  mechanisms.  There  has  been  evidence  for  the  predominance 
of  diffusive  effects  at  the  shorter  times  and  oxidation  for  longer  aging  periods  with  complete 
oxidative  control  at  the  higher  temperatures  (References  15  and  18). 

An  encouraging  correlation  was  found  for  the  tensile  strength  at  45°  to  the  fabric  warp  for 
CTL-91LD/Glass  Fabric  (Figure  13).  The  empirical  result  was 

TS(45°)  =  c1  exp  (c2w)  (23) 

A  strength  retention  of  only  about  10  percent  of  the  average  room  temperature  control 
(25,260  psi)  existed  for  w  =  0.8.  The  empirical  coefficients  and  other  results  have  been 
summarized  elsewhere  for  CTL-91LD/Giass  Fabric  and  a  second  material  (Reference  19). 
Other  attempts  at  mechanical  property  correlations  were  not  satisfactory  if  the  strength 
increased  with  an  increase  in  either  temperature  of  time.  Correlations  may  have  been  possible 
over  restricted  temperature  or  time  regions  as  suggested  by  the  similarity  of  the  plots  of 
Figure  14. 

Other  isothermal  air  aging  kinetic  experiments  for  aCTL-91LD  polycondensate  powder  and 
CTL-91LD/Glass  Fabric  were  informative.  A  -325  sieve,  100  mgm  powder  sample  had  a 
life-time  of  two  hours  at  600°F.  The  life-time  for  a  -325  sieve  100  mgm  laminate  powder 
sample  was  48  hours,  implying  a  glass  shielding  effect.  A  100  hour  life-time  of  a  1/8  inch 
thick  laminate  plate  was  undoubtably  associated  with  surface  area  and/or  volume  influences. 
These  results  suggested  a  need  for  a  flexible  empirical  kinetic  model  properly  accounting 
for  surface  area,  thickness,  two  controlling  mechanisms  with  a  transition  region,  and  probably 
catalysis,  oxygen  concentration,  and  weight.  An  established  model  would  be  necessary  prior 
to  any  serious  attempt  to  use  this  approach  for  the  analysis  and  correlation  of  the  mechanical 
properties  of  reinforced  laminate  plates  during  long-term  air  aging. 
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SECTION  VI 

SUGGESTIONS  FOR  FURTHER  WORK 

Additional  work  seems  necessary  to  establish  ranges  of  environmental  validity  for  the 
empirical  kinetic  model.  Particularly  questionable  aspects  include  heating  rate  extremes, 
isothermal  thermogravimetry,  micro  size  powder  and  large  laminate  plate  samples,  purging 
flow  through  a  powder  sample,  and  reduced  and  high  pressure  (Reference  27).  Current 
applications  for  hyperthermal  ablation  computer  programs,  mechanical  property  correlations, 
vacuum  life-time  predictions  make  realistic  environmental  simulation  an  essential  need.  Table 
XIII  summarized  some  conditions  pertinent  to  the  two  above  requirements  yet  within  practical 
experimental  capabilities  in  the  near  future.  Improvements  in  experimental  accuracy  of  up 
to  an  order  of  magnitude  better  than  for  the  accuracy  values  of  Table  X  were  desirable 
although  by  and  large  beyond  the  immediate  state  of  the  art  of  thermogravimetry. 

Improved  computer  programs  appeared  desirable  prior  to  further  study.  One  outstanding 
approach  gave  best  empirical  kinetic  parameters  in  a  least  mean  squares  sense  (Reference  43). 
Additional  empirical  kinetic  models  appeared  necessary  for  some  critical  cases  as  an  oxidative 
environment  or  sample  other  than  a  fine  polymer  powder  (Table  XIV).  Particularly  desirable 
features  for  future  computer  programs  included:  (a)  constant  heating  rate  and  isothermal 
versions,  (b)  evaluation  of  n,  A,  E  in  a  least  mean  squares  sense,  (c)  evaluation  of  confidence 
limits,  variances,  and  significant  digits,  (d)  experimental  corrections,  (e)  statistical  data  point 
weighting  for  the  entire  thermogram,  and  (f)  temperature  variance  options  for  n.  A,  E,  either 
singularly  or  in  any  combination. 

There  were  a  variety  of  polymers  for  which  definitive  kinetic  data  could  prove  helpful  in 
ablative,  coating,  structural,  and  other  aerospace  uses.  Polymer  classes  of  special  interest 
were  summarized  by  Figure  15.  A  future  goal  was  kinetic  information  on  the  base  polymers  as 
well  as  materials  using  such  reinforcing  agents  as  fibrous  carbon,  glass,  graphite,  and  silica. 
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APPENDIX  I 
GENERAL  THEORY 
The  constant  heating  rate  kinetic  relation 

-w  =  (  wnA/T)  exp  (-X ) 

was  rearranged  as 

/w  _n  .  rT  — y 

w  <iw  =  (A/T)  J  e  dT 

l  0 

for  the  assumption  of  a  zero  rate  at  zero  degrees.  K,  a  convenient  index,  became 

K  =  (  I  -  w )  for  n  =  0, 

=  -  In  w  =  l , 

=  (w1-n-l)/(n-l)  *  0,1 

Integration  of  the  other  side  of  Equation  1-2  by  parts  in  two  different  ways  gave 
K  =  (AE/RT)  p(X) 

p(  X)  =  X_le"X-/  X_le“XdX 
X 

.  r  •  _  m  +  l  -m  -2  -X 

K  =  (AE/RT )[  I  -  £  (“I  )  X  (m  +  I ) !  J  X  e 
1  m 

For  m  =  o,  a  relative  error  r  v/as  defined  as 


Therefore 


r  =  X2eX  p( X) 


K  =  .r  ( AE/RT)X  %  X 


Figure  9  presented  K  =  f(n,  w).  Figure  10  depicted  p(X),  r  =  f(X). 
The  maximum  rate  occurred  when 

d(-w)/dt  =  d(~w)/dT  =  d(~w)/dw  =  0 
Differentiation  of  Equation  1-1  with  respect  to  time  gave 


[nww°  *  +  (77T)wnx]  =  0 

L 


The  maximum  rate  followed  as 


'WTm  =  Wm(x/nT)m 


(I-l) 

(1-2) 

(1-3) 

(1-4) 

(1-5) 

(1-6) 

(1-7) 

(1-8) 

(1-9) 

(1-10) 


(I -ID 


11-12) 

(I- 13) 
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Also 

km  =  w^"n  (X/nT)mf  =  Ae~  (1-14) 

where  n  cannot  be  zero. 

The  two  principal  relations  for  the  maximum  rate  method  came  from  rearrangements  of 
Equations  1-13  and  1-14 

E  =  -W-  nR(TZ/w)  (1-15) 

Tm  m 

I  “*  (1  •  Xrr\ 

A  =  wm  (TX/nT)  em  (1-16) 

m  m 

The  third  necessary  relationship  for  n  resulted  from  a  comparison  of  Equations  1-10  and  1-14 


For  n  =  1 


and  for  n  /  0,  1 


l-n  , 

K  =  /n 

m  m 


=  e 
m 


w  =  f  r  {  n  '  -  I )  +  I )  1 
m  l  ■* 


l/(  n-l) 


{ I  — 17 ) 

(I- 18) 

1-19) 


It  should  be  noted  that  Equation  39  of  Reference  17  is  in  errQr  and  should  be  identical  to 
Equation  1-19  above. 
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APPENDIX  n 

COMPUTER  PROGRAMS  OUTLINES 


I.  Maximum  Rate  Experimental  Method  (MAXRAX),  See  Figure  1. 

A.  Read  input:  R,  T  ,  T ,  T  ,  w  ,  w  -\v_  »  w 

*  m  m  m  oi  Tm  r 

B.  Read  tables: 

'•  '|,k  *  MV 

2-  *,  *  M"i  ''i1 II. 

C.  Read  integral  values  of  n.,  r.  for  w  from  table. 

&  1  i  m 


D.  Calculate  n. 

J 


In  f  r.  .  (  n.  .  —  I  )  +  !  1 
L  '.J  i.J  J 


n. 


In  w, 


+  I 


m 


(niiiik^.n 


E.  Calculate  E,  X 


m 


in  E/RTm 


F.  Read  r.  from  table. 

G.  Calculate  n^  from  Equation  II— 3. 

H.  Calculate  E  from  Equation  II-5.  Calculate  A 


A  =  w'"nk(TX/nkT)m  exp  (X  ) 


m 


k'  'm 


m* 


I,  Write  output: 

1.  Run  identification. 

2.  n,  ,  A,  E,  r.,  X  ,  T 

k  j  m  m 

J.  Construction  subroutine 

II.  Transfer  Resolution  Integral  Method  (TRIM).  See  Figure  2. 
A.  Read  input: 

n»r»R,Tm*  '  wm’  woi  *  -^Tm,wr 


(I-l) 

(JE-2) 


(1-3) 

(H-4) 

(n-5) 

(X-6) 


(3T-7) 


x7 
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B.  Read  table: 

wQ  =  f,  (8)  with  9p<  8  <  9q  . 

C.  Calculate  T 

T  =  9  +  273.16 

D.  Calculate  w 

w  =  ( w0  -  wr )/{  w0 j  -  wr ) 

E.  Calculate  K 

B  K  =  (  I  -  w  )  for  n  =  0, 

=  -  In  w  =  I , 

=  (w'“n  -  l)/(  n  -I  )  *0,1 

F.  Calculate  y,  z 

y  =  log10(K/T2) 

z  =  l/T 

G.  Calculate  A,  E  by  linear  pseudo-method  of  least  squares 

y  =  bz  +  a 

H.  Calculate  A' 

A‘  s  ■'iTmw"ntm  e*P{Xfn) 

I.  Write  output: 

1.  Run  Identification. 

2.  n,  A,  E,  r,  T,  A’ 

J.  Construction  subroutine. 

III.  Parabola  Resolution  Integral  Method  (PRIM).  See  Figure  3. 

A.  Read  input:  n,  r,  R,  Tm,  T,  Tm,  wm,  wq.,  -wrfm,  wr 

B.  Read  tables: 

1.  wQ  =  ij (9)  with  6p<  8  <  8^. 

2.  nn ,  n0  .  .  .  n.,  with  i  =  6  maximum. 

C.  TRIM  subroutine. 


(3-8) 

(n-9) 

(I- 10) 

(3E-II) 

(1-12) 

(3-13) 

(3-14) 
(3  -15) 

( 3- 16) 

(3-17) 


(3-18) 

(3-19) 
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D.  Calculate  g 

g  =  I(woc  -  V2 

E.  If  g  =  0,  n  =  n' .  Proceed  to  item  (I). 

F.  Calculate  aQ,  bQ,  cq  by  parabolic  pseudo-method  of  least  squares 

Q  =  a0n2  +  bQn  +  cc 

G.  If  a  <  o,  designate.  Proceed  to  item  (V). 

H.  Calculate  n' 

n'  =  -  b  /  2  <r 
o  o 

I.  TRIM  subroutine. 

J.  Calculate  g ' 

9  =  I  <woc  ~  wQ)2 

K.  Calculate  w~ 

o 

*  [l.oj/q 

o 

L.  Calculate  s  (w  ) 

6  O 

se2(wo}  s  -2) 

0 

M.  Calculate  s  (w  ) 

pv  o' 

8p(V  s  [Z  (w0  -  W^)2] /q  -  |) 

N.  Calculate  F(wq) 

F(w0  )  =  Sp2  (w0)/s42(w0) 

O.  Calculate  a,  b 

a  =  log J0  r  ( AR/ET) 
b  =  -E/(  R  In  10 ) 

P.  Calculate  1/T 

T7T  =  [T  I / T j  /q 

2 

Q.  Calculate  s  (a) 

e  ' 

se(o)  = 


(£-20) 

(1-22)  | 
1-23) 

(31-24) 

or- 25) 

(It-26) 

(31-27) 

(31-28) 

(3E-29) 

(31-30) 

(31-31) 

(31-32) 

(31-33) 
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R.  Calculate  sg  (b)  _ 

s2(b)  =  s*(w  )/  IM/T-  l/T) 

S.  Construction  subroutine. 

T.  Substitute  s  temperature  span  for  q  temperature  span  with 

0,  <  8  <  9S 

U.  Repeat  Items  a)  through  (S)  for  s  temperature  span. 

V.  Write  output: 

1.  Run  identification. 

2.  For  the  q  temperature  span:  n\  A,  E.^R,  T,  A  ,  g  ,  wQ» 
sj>0),  8®<w0).  F(wo),  a.  b,  l/T,  TW.  Sg(a).  »e(b); 

*oc  'a1®1 


(n-34) 


m-35) 


(E-36) 


3.  Item  (V.2)  for  the  s  temperature  span. 
IV,  Construction  Subroutine.  See  Figure  4. 

A.  Read  input:  n.  A,  E,  r,  R,  T.  Tm>  wQi,  wr 


R. 

Read  table:  p(X)  =  f^(X) 

(It -37) 

C. 

Calculate  K 

K  =  ( AE/RT)  p  ( X) 

(H-38) 

D. 

Calculate  w 

w  =  1  “  K 

for  n  =  0, 

(3L-39) 

=  exp ( - K ) 

.  I/C  1  — n) 

=  [x(n-l)  +  !  J 

s  •  , 

*  0,1 

(It -40) 

(I -41) 

E. 

Calculate  wQC 

*oc  =  w^w0j  -wr)  +  wr 

(  1-42) 

F. 

Calculate  8 

B  =  T  -  273.16 

( 31-43 ) 

G. 

Write  output:  wQC  =  ^8) with  9r<  9  <  9g. 

(IE- 44) 
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V.  General  Notes. 

A.  Numerical  values  for  this  report  were: 

1.  j  =  k  =  30° 

2.  £  =  S  =  0.001. 

3.  w0,-l. 

4.  8,8  were  for  the  j  increment  nearest  to  w  ~  0.95,  w  ~  0.05, 

p  q  Pi 

respectively.  See  Equations  14-16. 

5.  0  =  140°,  6  =  1010°. 

r  s 

B.  The  following  was  omitted  for  brevity: 

1.  Standard  tests  for  avoiding  In  0,  In  (-number),  etc. 

2.  For  MAXRAX, 

a.  Subscript  m  omitted  for  n. 

b.  Subscripts  m,  i,  j,  k,  omitted  for  A,  E,  X  . 

m 

c.  A  manual  test  assuming  n  =  1  if  n^  k  <  0. 

d.  A  test  to  insure  that 

was  positive. 

e.  Rules  for  choosing  +  e  or  -  * . 

f.  -w„  ,  w  calculated  manually. 

im  m 

3.  Linear  least  mean  squares  pseudo-method  subroutine  for 
TRIM  with  output  A,  E. 

4.  For  PRIM, 

a.  Parabolic  least  mean  squares  pseudo-method  subroutine  with 
output  n,  g,  gc,  aQ,  bQ,  cq,  ±(gc  -  g)m,  Sg(g),  Sp(g). 

b.  Output:  n,  A,  E  and  wq  =  f(0)  as  an  option  for  each  n  value. 

c.  Equations  11-26  through  II-2y  for  the  flow  diagram. 

C.  For  TRIM: 

1.  n  =  n  (MAXRAX),  n  =  n.  (PRIM),  or  n  =  n'  (PRIM). 

m  i 

2.  r  =  r  (MAXRAX) ,  and  T  =  T  unless  noted. 

m  '  m 


nr -45) 

(2-46) 

(2-47) 

(I-4t 

(H-49) 


(2-50) 

(2-51) 


(2-52) 
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D.  For  Construction  Subroutine: 

1.  Run  identification  and  n.  A,  E,  r,  X  ,  T  included  in  output 
when  run  separately. 

2.  s  temperature  span  except  as  a  special  case  for  PRIM. 

E.  Symbols  f  and  symbols  using  subscripts  i,  j,  k  usually  unique 
for  each  program  outline. 
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*Constant  rate  and  isothermal  mass  spectormetric  thermal  analysis  results  included  for  completeness  (Reference  51). 
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TABLE  VII 

PRIM  EMPIRICAL  KINETIC  PARAMETERS 
FOR  PURE  AND  GLASS  FABRIC  REINFORCED 
PHENOL-FORMALDEHYDE  POLYCONDENSATE  POWDERS 


A.  Sample  and  Method 


Atmosphere 

He  purge 

He  purge 

Code 

A-4  Casting 

A-2  Laminate 

Crucible 

Coors  230  3/0 

Coors  230  3/0 

Fabric  Reinforcement 

none 

yesa 

Grinding  Method 

cutter;  liquid  N2 

hole  saw;  liquid 

Particle  Size 

-325  sieve 

-325  sieve 

Postcure  Temperature,  b°C 

178 

178 

Resin  Content,  % 

100 

27.0 

Trade  Name 

CTL-91LD 

CTL-91LD 

Weight,  mgm 

20 

100 

Empirical  Kinetic  Parameters 

n 

.1.62580 

1.73258 

A,  min  _1 

24.234 

201.572 

E,  kcal/mole 

9.18206 

12.9362 

q  span,  °C 

350-770 

350-770 

rm 

0.787766 

0.845028 

1,  °C/min 

7.7 

7.7 

Tr,  °C 

950 

950 

wr 

0.552 

0.8888 

C.  Conformance  Indices 


g'  x  103 

s2(a)  x  105 
e 

s2(b)  x. 102 


0,5973 

0,2667 

0.2347 

0.009806(620°) 


0.05057 

0.02258 

-0.01987 

0.003142(620°) 


a112-"b''  glass  fabric  with  A1100  finish. 
^Maximum  postcure  temperature. 
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°  181-"E”  glass  fabric  with  A1100  finish. 
b  About  1/2 x  1/4 x  1/8  inch. 
c  Estimated  maximum  postcure  temperature. 
d  Estimated  variation:  25.85  -  3.25 

•Method  of  Freeman  and  Carroll.  See  Table  II  for  results  for  an  empirical  method. 
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TABLE  IX 

MISCELLANEOUS  PRIM  RESULTS  FOR 
PHENOL-FORMALDEHYDE  POLYCONDENSATE  POWDERS 


A.  Sample  and  Method 


Atmosphere 

Ar  purge 

Ar,  He  purge 

Code 

B-C 

See  Trade  Names 

Crucible 

porcelain 

Coors  230  3/0 

Grinding  Method 

Liquid  Nitrogen 

cutter,  drill 

Room  Temperature 

filing 

compaction,  filing 

Particle  Size 

fine 

fine,  -325  sieve 

Postcure  Temperatures,  °C 

230 

178,204,230 

Trade  Name  (Or  Code) 

B-C 

B-C,  CTL-91LD,  RI4009,  SC1008 

Weight,  mgm 

200 

20,  100,  200 

Empirical  Kinetic  Parameters 

n 

1.61375b 

1.55578+0.07002 

A,  min 

51.5221 

60.9338+41.3902 

E,  kcal/mole 

11.4162 

10.38033+1.19827 

q  span,  °C 

350-830 

560  ±300 

Yra 

0.780755 

0.7924905  ±0.0117355 

T,  °C/min 

3.3 

5. 125 ±2. 575 

Tr>  °C 

815 

882168 

wr 

0.555 

0.577+0.025 

C.  Conformance  Indices 

g'  x  103  0.6770  0.37462+0.30238 

s2(a)  x  10s  0.2351  0.14589+0.12081 

s2(b)  x  102  -0.2115  -(0.128635+0.106065) 

±(w  -  w  )  0.01694(440°)  0.011148+0.005792 

oc  o  m 

D.  Reference  Lochte,  et  al  (37)  This  Reports 

Maximum  processing  temperature. 

^Approximate  estimates  for  q  span,  r_,  vr.  w0  values  from  a  smooth  curve  from  tabulated 
data  extrapolated  graphically  to  950°.  A'  was  not  calculated  accurately. 

Averages  from  Tables  VI  (3°,  6°,  9°/min),  VII,  VIII  and  IX. 
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SOME  ADDITIONAL  LITERATURE  EMPIRICAL  KINETIC  PARAMETERS  FOR  PHENOL- FORMALDEHYDE  POLYCONDENSATE  MATERIALS 


Ai,ML-TR-65-24G 
Pirt  II 


DO 

c 

—  1 
<r>  & 

S.U.  o 


^  w 
03  *2  m 

^OQ- 


>  to 

*  CD  CO 


•  < 

X 

'fe 

5s* 


<D  03  «S 
O  ■—  <13 

q_lu.ce: 


a 

eg 

DO 


-A  ^ 

•— «  O  CD 

cd 


2  S 

20- 


X 

CD 

O^rco 


—  5 

ro  .£ 


o 

x 

03 


JZ  CO  CL. 


o 

o 


X3 

03 

g.s 

,p  v> 

:.£<g 

:S« 

~  </>  2 
P  co  £ 

J2_g 

CJ  Q_ 


I 


<a 

y 

"a. 


22  ~ 


03 

■S 

o 


s 

1  i&s 

iz  o  . 

2  S-S 

if  =£2 

-  '« g 

?  “<5 

o  _ 

CL  'C-— 
c=  _S  t/i 


Q_ 

o 

03 

C 

SZ 

15 

o 


o 

_  B 


B  .2 

cs<  Ui 


-  c/3  — 

03  <V 
DO  O  c= 


£ 


<g  q: 
f= 


O  "O 

-  ft 

</)  =3 

ft  <* 

15  J 

>  <g 
•o  > 

DO  LCl 


< 

_D 


41 


SOME  PROBATIVE  KINETIC  EXPERIMENTAL  RANGES  FOR 
POLYMERS  AND  REINFORCED  PLASTICS3 «f 


AFML-TR-65-246 
Part  II 


a 

B 

B 

B 

00 

oo 

60 

60 

• 

Li 

a 

B 

a 

B 

B 

0 

43 

B 

S 

00 

B 

0 

6 

u 

60 

o 

o 

60 

ob 

a 

60 

o 

60 

o 

CO 

tH 

o 

o 

o 

o 

rH 

0! 

CM 

04 

in 

m 

ITI 

m 

04 

m 

04 

m 

1 

3: 

1 

| 

04 

04 

1 

04 

1 

04 

1 

04 

m 

H 

H 

*H 

H 

r4 

H 

H 

rH 

rH 

• 

• 

i 

i 

• 

I 

• 

1 

• 

1 

o 

o 

o 

H 

* — 4 

O 

H 

o 

r-4 

o 

*H 

• 

o 


i 

>4 

0)  C 

U  -H  c\l 


04  ^ 

M 

04 

04 

1 

04 

1 

H 

• 

rH 

o 

<D  <D 

04 

rH 

» — I 

*H 

04 

<>4 

04 

04 

O 

o 

o 

60  AJ 

| 

• 

• 

• 

1 

1 

1 

1 

1 

i 

vO 

^  r! 

3 

0< 

rH 

O 

O 

o 

H 

r4 

H 

H 

o 

o 

1 

m 

«H 

W 


CD 

#H 

H 

rH 

rH 

rH 

u 

1 

i 

1 

1 

1 

oo 

CO 

• 

3T3 

• - I 

rH 

*H 

• — i 

tH 

»H 

rH 

rH 

1 

1 

m  a 

o 

1 

i 

i 

o 

o 

o 

o 

OI 

o\ 

X 

W  AJ 

o 

m 

in 

m 

o 

o 

o 

o 

1 

1 

• 

cu  ca 

• 

• 

• 

• 

• 

• 

• 

o 

o 

AJ 

#k  * 

o 

o 

o 

o 

o 

o 

o 

o 

H 

rH 

3 

to 

IX 

o 

I 

d 

M 


-  O 
fX  4> 


3 

x> 


-  <u 

N  4-1 


rH 

a) 

O 

o 

o 

d 

d  d 

CO 

u 

o 

O 

o 

m 

o 

o 

O 

o 

•A 

-H  -H 

E 

0 

o 

o 

o 

CO 

vO 

vO 

O 

o 

in 

HH 

p 

aj  a 

vO 

vO 

vO 

rH 

t — i 

H 

vO 

vO 

aJ 

d 

<u 

CU  O 

1 

1 

1 

l 

1 

1 

I 

I 

M 

O  43 

.c 

M  o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

~A  U 

AJ 

(1) 

o 

o 

o 

m 

m 

m 

o 

o 

•a  <y 

U  CO 

o 

Or 

<r 

<r 

•H 

«H 

rH 

<r 

<r 

<u  > 

• 

cO  O 

co 

B 

Li  -A 

CO 

iH  U 

M 

CD 

O  X> 

CD 

M  d. 

H 

(3  •  a) 

iH 

CO  (X 

o  -a 

CO 

MO  -rl 

AJ 

M  O  X 

0) 

O  O 

a)  in  o 

EH  U 

60 

,o 

rH  cn 

o 

4J 

C 

*  C 

o 

C  — 1  »H 

CD 

cu 

*H 

0)  -H 

o 

O  o 

oo  a  oo 

4J 

aj  E 

VO 

vO 

vO 

vO 

vO 

vO 

vO 

vO 

o 

CD  *H 

CO 

(U  \ 

1 

i 

I 

1 

1 

1 

I 

I 

U  Li  c 

M 

AJ  X 

<D 

0$  O 

CO 

CO 

C0 

rH 

rH 

rH 

CO 

CO 

a)  o 

CD 

CO 

3: 

o 

•a  <u  -h 

XJ 

Sn 

H  >  U 

rH 

CO  (D 

O  iH  n! 

o 

AJ 

43  Li  N 

42 

M  cu 

(0  iH 

O  H 

& 

a)  13  rH 

CD 

«A4  Dh 

u 

u 

u 

u 

rH  -fH  T-i 

rH 

4J 

0) 

<d 

0) 

a) 

<D 

(D 

CD 

a> 

CD 

a) 

CL  X  43 

CU  X>  <D 

a) 

XJ 

T3 

AJ 

AJ 

4J 

T3 

AJ 

XJ 

4-4 

EOrt 

• 

a 

Q)  44 

8 

3 

cU 

CU 

S 

cU 

3 

cU 

3 

(U 

3  44 

B 

to 

d  -rl 

o 

o 

O 

rH 

fH 

O 

rH 

O 

rH 

O 

rH 

M  •  M 

AJ 

M 

oo  d 

a> 

04 

0- 

&4 

04 

04 

04 

04 

04 

04 

04 

U 

CU 

•rl  -H 

o 

Co  L) 

d 

M  OH 

<U  O  43 

rH 

0) 

0)  d 

OO  600 

0*0  H 

OiOtI 

o 

o 

1 

rH  (1) 

* 

M  -rl 

T3  so 

M 

w  a 

o 

c  o 

3 

d  3 

sun 

AJ 

CO 

0)  M 

CD 

a) 

o 

CU 

M 

a 

<d 

*4 

rC 

<u  a)  iw 

QJ 

•rl  43 

0) 

a 

M  JH 

<w 

> 

M  O 

rO 

CO 

U 

a)  a)  <u 

AJ 

<0  co 

cU 

Or 

o 

4J 

43  43  e 

cu 

T3 

D.  O 

•H 

CO 

B 

a) 

(X  CL  tH 

u 

0) 

X  d 

M 

O 

aj 

<D 

e 

in  m  l) 

d 

0)  cx 

(U 

B 

< 

u 

-H 

o  o 

0) 

•H 

cx 

> 

aj 

rH 

> 

6  3  fn 
u  u  p 

60  4H 

Q)  co 

< 

<D 

o< 

cu 

cU 

0) 

u 

c 

Li 

> 

6 

V4 

<D 

CO  (3  rt 

o 

(0  O 

O 

a 

•H 

60 

U 

60 

in 

cu 

o 

rH  Li 

*r-> 

aj 

•H 

AJ 

a 

0) 

O 

3 

Li  u  m 

CL 

CU 

rC 

B 

CU 

*H 

rC 

CO 

u  u  <D 

o 

o 

l! 

s 

60 

cd 

T3 

AJ 

AJ 

P 

CO 

CL)  (U  O 

u 

CO 

rA 

•H 

<2 

*H 

cU 

o 

<D 

(U 

G  G  <D 

0) 

» — i 

rH  S 

<D 

£ 

X 

<D 

CO 

J-» 

H  H  21 

N  <  <  o 

3: 

o 

H 

H 

04 

<U  Xi  O  TU  (D  MH  »H 

42 


SOME  PROBATIVE  KINETIC  RELATIONS  FOR  EMPIRICAL  INVESTIGATION 
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Figure  1.  Flow  Diagram  for  the  Maximum  Rate  Experimental  Method 
(MAXRAX)  Computer  Program 
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Figure  2.  Flow  Diagram  for  the  Transfer  Resolution  Integral  Method  (TRIM) 
Computer  Program 
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Figure  5.  Thermograms  for  Phenol-Formaldehyde  Molding  Powder 
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Figure  6.  PRIM  Parabolic  Error  Plots  for  Phenol-Formaldehyde 
Molding  Powder 
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Figure  7.  Thermograms  for  Phenol-Formaldehyde  Polycondensate  Samples 
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Figure  12.  Reduced  and  Calculated  Thermograms  for  CTL-91LD/Glass  Fabric 


PERCENT  OF  ROOM  TEMPERATURE  STRENGTH 


RESIDUAL  VALUE  -  (PERCENT  OF  ROOM  TEMPERATURE  CONTROL) 


AFML-TR-65-246 
Part  II 


TEMPERATURE,  *F 


Figure  14.  Strength  and  Weight  Retention  for  100-Hour  Air  Aging  of  CTL-91LD/Glass 
Fabric 
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